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The 'Ghost' interference experiment is analyzed when the source of entangled photons is a mul- 
timode Optical parametric Amplifier(OPA) whose weak limit is the two-photon Spontaneous Para- 
metric Downconversion(SPDC) beam. The visibility of the double-slit pattern is calculated, taking 
the finite coincidence time window of the photon counting detectors into account. It is found that 
the coincidence window and the bandwidth of light reaching the detectors play a crucial role in the 
loss of visibility on coincidence detection, not only in the 'Ghost' interference experiment but in 
all experiments involving coincidence detection. The differences between the loss of visibility with 
two-mode and multimode OPA sources is also discussed. 



PACS numbers: 

I. INTRODUCTION 



The 'Ghost' interference experiment is typical of two- 
photon interference experiments that bring out quantum 
entanglement features of light. These features are not 
just confined to the appearance of the double-slit pattern 
on coincidence detection. The independence of the result 
on where the optical elements are situated in the experi- 
mental set-up (a double-slit in the case of 'Ghost' inter- 
ference) is a key feature of entanglement. The 'Ghost' 
interference experiment has been performed in the low- 
gain limit of parametric down-conversion 0. It is of in- 
terest to examine this experiment in the high gain regime 
of parametric amplification. Analysis of similar experi- 
ments using two-mode OPA states 2| suggest the loss 
of visibility at large gains. Recently different detection 
schemes have been proposed to circumvent this prob- 
lem 0]. Here we present a detailed calculation of 
the 'Ghost' interference experiment using a multimode 
OPA as the source of entangled photons. We analyze 
the effect of the coincidence time window of the photon 
counting detectors on the experimentally observable in- 
terference pattern and visibility. We attempt to explain 
how the properties of the source of entangled photons and 
experimental limitations affect the observable interfero- 
metric effects. We find that the loss of visibility depends 
on the coincidence time window and the bandwidth of 
the source. The coincidence time window causes a loss in 
visibility at much lower gains than expected with ideal 
detectors. A multimode source causes loss in visibility 
even if ideal detectors are used. Combined with a fi- 
nite coincidence time window, a multimode source can 
reduce visibility to less than 0.5 even at very low para- 
metric gain. These effects are not specific to the 'Ghost' 
interference experiment but occur in all experiments with 
similar sources and detection schemes. 
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II. MULTIMODE INTERACTION IN OPA 

We consider a non-degenerate OPA Q comprising a 
non-centrosymmctric crystal with a second order non- 
linear susceptibility ■ The crystal is pumped by a 
CW laser, given by 

E p (r,t) = EM+El-t (1) 
= E p e l ^- p -^ t] + H.c (2) 

where H.c stands for Hermitian conjugate. We assume 
that the pump is not depleted, has a constant amplitude 
E p and can be treated classically. If a weak multimode 
signal is also input into the OPA, the output after para- 
metric interaction within the crystal consists of amplified 
modes of the signal and idler where the idler modes obey 
the phase matching conditions, 

LU p = LU S + UJi (3) 

hp — k s + hi (4) 

where k p , k s and ki are the wave vectors of the pump, 
signal and idler photons. The quantized field of the signal 
and idler inside the crystal are given by Q , 

= T.^^-^H.+H.c (5) 

kg 

k s 

Here is the annihilation operator for a photon in mode 
k. Note that the dispersion relation oj(k) is different from 
that in vacuum, u>(k) ^ ck. 

The non-linear interaction between the crystal, the 
pump, signal and idler modes are characterized by the 
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interaction Hamiltonian |j| , 
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where the integration is over the volume of the crystal 
illuminated by the pump. For a long crystal and a wide 
pump beam, 
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The interaction Hamiltonian is simplified to 
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The delta function in Eq. JJjJ indicates the entanglement 
of the OPA states in wave vector space. Combining the 
photon commutation relations and the time evolution 
equations for signal and idler modes, the time-evolved 
signal and idler are given by |ld , 

= a kS°) cosh (\tk s \) ~ ial (0)sinh(\Z ks \) (10) 

a\{t) = at(0)cosh(^ k j) + ia^(0)sinh(\^ ks \)(ll) 

where t is the average time taken by the photons to cross 
the crystal, (0) and at (0) are the annihilation and 
creation operators for the input signal and idler mode. 
The factors cosh(\£k s \) an d sinh(\£k,\) are the amplifi- 
cation factors that depend on the strength of the non- 
linearity %( 2 ) ,the pump E p and the frequency of the sig- 
nal and idler modes. Perfect phase-matching ensures that 
each signal mode interacts with only one idler mode, se- 
lected by the phase-matching conditions. Eq's. IjlOl) and 
(|ll|l relate the photon creation annihilation operators at 
the output face of the crystal to those at the input. They 
are derived from the unitary transformation, 
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Here £h s — (x^ 2 ' ] E p ^uj s ujit)e™/ 2 , may be recognized as 
the squeeze parameter of the OPA . The state at the 
output of the crystal is given by 
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(13) 



which is a multiphoton state. In the limit <C 1, the 
expansion of Eq. (|13fl can be limited to first order in the 
interaction Hamiltonian, giving a vacuum state and an 
entangled two-photon state. This is the two-photon limit 
of SPDC. 



III. EXPERIMENTAL SET-UP AND 
CALCULATION 



— ► 




OPA 






I 




/ 








PBS 













FIG. 1: Schematic of the experimental setup for observing the 
two-photon 'Ghost' interference. For simplicity, the prism to 
remove the pump is not shown. The double-slit has width 
b = 0.165 mm and slit distance a = 0.4 mm. The relevant 
distances in this experiment are d\ = 0.3 m, = 1.5 m, and 
d 1 = 1 m. 



We now look at the experimental set-up of the ghost 
interference experiment Q. A non-linear crystal is 
pumped by a CW laser(A p = 351nm), to generate pairs 
of collinear, orthogonally polarized signal(e-ray) and 
idler(o-ray) photons (Type-II SPDC). The signal and 
idler beams are separated by a polarizing beam splitter. 
The signal beam passes through a double slit aperture 
to a photon counting detector D\. D\ is a fixed point 
detector. The idler is scanned by an optical fiber and 
the output from the fiber is coupled to a detector (D 2 ). 
fi and ji are filters (spatial or spectral) that limit the 
number of wavevectors (and bandwidth) of light reaching 
the detectors. The detectors in the signal and idler arms 
are connected to a coincidence circuit. 



IV. COINCIDENCE DETECTION 

We are interested in coincident detection from the sig- 
nal and idler when the input state is a vacuum. The 
probability of detecting a photon in D\ in position r\ at 
time t\ and another in D 2 at r*2 and t 2 is proportional to 
the second order correlation function [jj , given by, 

G^\r l ,t 1 ,r 2 ,t 2 ) = {^) (14) 



\*)=E { +\r 1 ,t 1 )E [ 2 + \r 2 M)\0) (15) 

where E\(rx,ti) and E 2 (f 2 ,t 2 ) are the electric fields at 
the detectors D\ and D 2 . The fields at the detectors can 
be expressed in terms of the fields at the output face of 
the crystal through Green's functions which describe the 
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propagation of the beams through the optical system. 
The positive frequency part of the fields at the detectors 
are then given by 



k 

E ( 2 + \r 2 ,t 2 ) = 
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Cl gA(k ,n)a ko (16) 



2 gB(k e ,f 2 )a kc (17) 



We now look at the state of the fields at the two detec- 
tors when the input state is a vacuum. Note that we are 
relating the fields at the detectors to the vacuum at the 
input face of the crystal. This involves two transforma- 
tions - from the detectors to the output face of the crystal 
through the Green's functions and from the output of the 
crystal to its input by the unitary transformation of the 
operators in Eq's. (fTUf) and ifTTJl. 

Using the operator commutation relations we get, 
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From a careful examination of the above equation, we 
find that the first term indicates correlation between the 
modes detected in the two arms. The second term has 
no such correlation since it factors into two independent 
terms, one for each arm. Physically, this term corre- 
sponds to accidental coincidences of photons that are not 
entangled with each other. We expect that the double slit 
diffraction pattern will emerge from the correlation in the 
first term while the second term causes a loss in visibility. 
The second order correlation function G^ can now be 
written as, 
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Implementing all the 6 functions and converting the sum- 



mations in G^ 2 ) to integrals > we have 
G< 2 > = 
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The expression for G^ 2 -* in Eq (|18fl can be rewritten as 
follows to emphasize the nature of the correlation in each 
term. 
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where G^ and G B ' are the first order correlation func- 
tions for the two arms, the The coincident counting rate 
is then calculated using, 
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To make the discussion easier, we shall evaluate the two 
terms in Eq 1)20(1 separately. First we make a few ap- 
proximations. It is usually easier to work with quantities 
invariant along the beam: w, the angular frequency, and 
q, the component of the wave vector parallel to the out- 
put face of the crystal. The z component of the wave 
vector for a photon of polarization (3 is 



K = 



(21) 



Outside the crystal, the index of refraction np = 1. Inside 
the crystal, np depends on the orientation of the optic 
axis with respect to the wave vector of the beam 
We assume \q\ <C |fc| and k z = ujnp(u)/c, so that 



dk z — — 

up 



(22) 



where up is the group velocity of a photon of polarization 
f3. Since the integrals in Eq. I|18f) are over modes outside 
the crystal, 

Jfk^jfgto^J*^ (23) 

We make the following approximations to further sim- 
plify the problem. We assume that the central frequen- 
cies of the signal and idler are degenerate, ie 

= Vl + v (24) 

= fi e - V (25) 
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Further we will assume that v <C f2 , e and uj = uu e = 
lu p /2 everywhere except the exponential terms. The 
squeeze parameter | is considered constant for all 
modes. With these approximations, we can now calculate 
the first term in the rate of coincidence counting. 
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Completing the time and frequency integrals in Eq. 132( 1 , 
we can now infer the effect of the coincidence window T 
in a general case without entering into the details of the 
experiment. 



Rent represents 'true' coincidences of mutually entan- 
gled photons. The Green's functions in Eq. 1(28(1 imply 
that that the detector Di is fixed at the origin so that 
f\ = (d\ + di)z and the position of detector D 2 is 
t 2 = p 2 p + d 2 z(see fig 1). T is the time window of coinci- 
dence detection. If the two detectors register two photons 
within time T of each other, the photons are assumed to 
be part of one entangled pair. For the present calcula- 
tion T = 1.8 ns. We have also used the fact that for a 
given pair of signal and idler modes related by the phase 
matching conditions, q Q = —q e . Since the visibility, being 
a ratio, is not affected by the finite detection area of the 
detectors, we will not bother about them in this calcu- 
lation. We now look at the frequency and time integrals 
in Eq. 1(27(1 . If the bandwidth of light reaching the two 
detectors, A, is such that AT 3> 1, we can approximate 
the frequency and time integrals by 2ttAT y| lljj ■ 



Rent = 2ttA cosh 2 (\Z\)sinh 2 (\Z\) (29) 

J d 2 qeg A (^Y, ~q e , 0, di + d\)g B (^-, q e , p 2 ,d 2 ) 



Using the Green's functions from the appendix A, R en t 
(ignoring common constants) is given by, 
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where t(q) is the fourier transform of the aperture func- 
tion in wave vector space. So for a double-slit aperture 
we have the expected 'Ghost'interference pattern in the 
true coincidences. cosh 2 (\^\)sinh 2 (\^\) is the amplifica- 
tion factor. Now we look at the accidental coincidence 
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Comparing R ent and R acc from Eq's. (|2"§|) and (J2SJ, we 
find that at low gain, R acc is negligible compared to Rent 
since smft, 4 (|£|) — > (|£|) 4 - The coincidence window has 
no effect on the result as only a single pair of photons is 
produced within the time T. But as the gain increases, 
and more pairs of photons are produced within a given 
time interval, Race becomes significant and the visibility 
of the interference pattern begins to fall. 

Using the Green's functions and a double-slit with 
aperture function, 
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where a and b are the distance between the centers of 
the slit and width of the slits respectively and B is the 
length of the slit, R aC c (ignoring common constants) for 
the 'Ghost'interference experiment is found to be 
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For a given gain value |£| and coincidence window T, 
Race is a constant since the interference pattern behind 
the double-slit in arm A (see Fig. 1) is averaged over due 
to the bandwidth of the source and the intensity at the 
scanning detector, uniformly illuminated by the SPDC 
beam, is a constant. The details of this calculation are 
given in appendix B. 

The complete expression of the coincidence counting 
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rate in the 'Ghost' interference pattern is given by 



R oc R ent + R ac 
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The visibility of the interference pattern as a function of 
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FIG. 2: The visibility of interference fringes as a function of 
parametric gain, when the coincidence time window T — 1.8ns 
and bandwidth of the light reaching the detectors is (a) lnm 
and (b)10nm 

the parametric gain, calculated from the expression for 
coincidence count rate in Eq. I|36[l . is given by 
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where A and B are constants arising from experimental 
factors. Fig. 2 shows plots of visibility as the bandwidth 
of light reaching the detectors is increased. 

We now look into the features of the visibility plots and 
analyze the factors that give rise to these features. The 
important parameters for this purpose are the gain terms 
cosh(\£\), sinh(\£\) and AT in Eq. (jHEJ- The bandwidth 
A is a measure of the number of modes since we have 
assumed perfect phase matching. AT can be thought 
of as a measure of the detectors' ability to resolve two 
entangled pairs of photons. A value of AT ^> 1 leads 
to increase in accidental counts since detectors cannot 
distinguish every entangled pair of photons. The effect 
of the coincidence window time T is easy to understand. 
In the very low gain limit (sinh\^\) w 0) only one single 
pair of entangled photons is produced within time T and 
so every entangled pair can be distinguished. As the gain 
increases, many more entangled pairs are produced and 
reach the detectors within the coincidence window time 
causing accidental coincidences of photons belonging to 
different entangled pairs, and hence visibility is lost. 

To understand the effect of the number of modes on 
the visibility, we go back to the output state of the OPA. 
The OPA state is restricted to include the two photon 
state and the next higher order interaction giving four 
photon states. 

The output at the OPA, given in Eq. (|T3)) can be writ- 
ten as 



n 



exp[ 



■4(0)oJ.(0)]|0) 



(39) 



where k s + ki = k p is the phase matching condition be- 
tween the signal, idler and pump modes. We consider n 
pairs of signal and idler modes, expand the exponential 
operator term in Eq. (|39|l and omit terms greater than 
second order in the gain parameter £. If the gain is con- 
sidered constant for all the modes, then the unnormalized 
OPA state is given by 



W = |o) 



(40) 
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From the expression for the truncated OPA state in 
Eq. (|40|l we infer that the states in the first and second 
terms of the equation lead to 'good' coincidence counts. 
The detectors detect photons belonging to a single en- 
tangled pair or to an entangled four photon state. The 
third term, on the other hand, can lead to detection of 
two photons belonging to different entangled states ('bad' 
counts) half of the time. The conditional probability of 
getting a good count, given n modes, is found to be 



P(good\n) 



16 + < 2 



(41) 



6 



For a large number of modes, this probability tends to 
0.5. This implies that when a large number of modes are 
allowed, the visibility of the coincidence detection pattern 
falls to 0.5. 

Though the number of good counts seem to dominate 
according to Eq. 141|) . the number of modes along with 
the coincidence window time produce a loss in visibil- 
ity greater than 0.5. As the gain increases the proba- 
bility of both good and bad counts increase and tend to- 
wards a constant limit. This leads to the flattening of the 
visibility with rising parametric gain. The flattening of 
the visibility occurs at lower gain as the bandwidth(and 
number of modes) increases. The value of the limiting 
visibility falls as the number of allowed modes increases. 
This is to be expected since as the gain £ rises, entangled 
pairs (the first order states of the OPA) are emitted in 
all modes. If the number of modes allowed in the experi- 
ment are increased then the number of possible good and 
bad counts also increase. The coincidence time window 
T further adds to the bad counts causing a greater fall 
in visibility and the limiting visibility is lowered. 

If a smaller number of modes are allowed into the ex- 
periment, for example by using a fine pinhole, the vis- 
ibility can be maintained at values greater than 0.5 for 
larger values of gain. But even in the case of just two 
modes, the visibility eventually falls due to the higher 
order states of the OPA and tends to a constant as 
cosh(\£\) « sinh{\^\) rj e^l/2. 

This is the fundamental difference in the mechanism of 
visibility loss with two-mode and multimode OPA. In a 
multimode OPA the loss of visibility is mainly due to the 
number of modes, and occurs at much lower gain than 
the two-mode OPA, where the higher order terms lead to 
loss of visibility for a given coincidence window time. 



V. DISCUSSION AND SUMMARY 

We have analyzed the effect of a multimode optical 
parametric amplifier source in the 'Ghost' interference 
experiment, taking into account the finite coincidence 
window of the photon counting detectors. We find that 
the loss of visibility with increasing parametric gain is 
strongly dependent on the coincidence time window. A 
longer coincidence time window reduces the ability of the 
photon counting detectors to resolve entangled pairs as 
the parametric gain of the OPA increases. We have also 
highlighted the differences between effects observed with 
a two-mode and a multimode source. An increase in the 
number of modes in the experiment increases the proba- 
bility of accidental coincidences between photons belong- 
ing to different entangled pairs (or states). Further, this 
experiment limited loss of visibility occurs in all coin- 
cidence counting measurements though it is significant 
only in the regime of a strong source of entangled pho- 
tons like an OPA. We conclude that a cautious choice of 
sources and detection schemes are necessary in order to 
observe certain signatures of entangled light in a macro- 



scopic regime. 

Since the completion of this calculation we have be- 
come aware of a two-photon absorption technique demon- 
strated by Dayan et. al. ^(| where Rb atoms undergoing 
simultaneous absorption of signal and idler photons over- 
come the problem of temporal resolution associated with 
a strong broadband source like the multimode OPA. The 
two-photon transition is sensitive to minute delays (order 
of 100 fs) between the signal and idler photons. But such 
a detection scheme does not discriminate between entan- 
gled and separable pairs of photons and cannot reduce 
the loss of visibility in a 'Ghost' interference experiment. 
Further, the use of such highly sensitive detection scheme 
requires a high gain OPA source which enhances acciden- 
tal coincidence counts. 
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APPENDIX A: GREEN'S FUNCTIONS FOR 
PROPAGATION THROUGH A LINEAR 
OPTICAL SYSTEM 

We give a brief review of propagation of an electric field 
through a diffraction limited linear optical system in each 
arm of the experimental set-up, following the treatment 
in 0] . The positive frequency part of the electric field 
at a time t at the input of a detector at r = zz + p is 
given by 

J dLud 2 qE(Lu)g(q,LU,p,z)a fj (LU,q) (Al) 

where ap(uj,q) is the annihilation operator at the source 
for a photon of angular frequency u>, transverse wave vec- 
tor q and polarization j3. The unit vector z, is the inward 
normal to the detector surface. E(u>) is a slowly varying 
function required for dimensional reasons and can be as- 
sumed constant in the current analysis. g(q, u>, p, z) is the 
optical transfer function or Green's function which de- 
scribes propagation through the linear optical system. In 
classical electromagnetic theory g connects electric fields 
in real space. In our quantum mechanical analysis, it 
connects two operators in photon number or Fock space. 
The superposition principles involved in calculating g are 
purely classical from classical electromagnetic theory. So 
in the quantum mechanical context it is best thought of 
as arising from boundary conditions on the modes of the 
fields irrespective of the state of the system. 

We now calculate the Green's function for the arm A in 
the Fig. (1). The Green's function is expressed in terms 
of the aperture function defined by t{p a ). 
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g A (q,u,0,di + d x ) 

j d 2 p a J d 2 p s h LU (-p a ,d' 1 )t(p a )h u] (p a - p s ,di)e^* 

(A2) 

where /o s and /9 a are the transverse co-ordinates of the 
source (crystal) plane and aperture plane. In the Fresnel 
approximation [17]. 
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The Green's function in the arm B of the experimen- 
tal set-up, for each plane wave mode, assuming that the 
source has a large cross-section is given by 



(A6) 



(q,uj,p,d 2 ) = d p s h u (p- p s ,d 2 )e lq '' 
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^(\q],—d 2 ) (A7) 
Using the above expressions for g A and g B , we see that, 

J d 2 qg A (~q,uj,0,d 1 + d[)g B {q,uj,p,d 2 ) = 

(^) e<( " /e)( y ' 1+d2) JdWGftl, ^(Pa) 

(A8) 

In the far field Fraunhofer approximation, the tp's in the 
above expression go to unity and 



In the far- field Fraunhofer approximation, 
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APPENDIX B: SINGLES DETECTION 



We now look at results of detection at the detectors 
D\ and D 2 individually, without caring for coincidences. 
This involves calculating the first order correlation func- 
tions 0, 



G^(f t ,U) = {Q\E[-\n,t t )E { +\n,u)\0) (Bl) 



where i — 1,2. Using the same techniques as in coinci- 
dence detection, we find that the first order correlation 
functions are given by 
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(A9) (r 2l t 2 ) is the correlation function in the arm without 



where t is the fourier transform of the aperture function. 
Further. 
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(A10) 



the double slit and as expected, there is no interference 
pattern at D 2 . The correlation function in the arm with 
the double-slit G^(fi,ti), suggests that there might be 
a pattern at D\ , if the transverse components of the wave 
vector (q a ) reaching the detector is narrow enough. For 
the wavelengths and the dimensions of double-slit chosen, 
an interference pattern will be observed if the detection 
angle of the SPDC beam at the detector A9± < 2X p /a w 
1.5 mrad. But for the lnm filter used, A8± « 15 mrad. 
Therefore, no interference pattern is found behind the 
double slit due to the large divergence of the SPDC beam. 
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